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ABSTRACT 

We report on the latest discovery of the HATNet project; a very hot giant planet orbiting a bright 
(V = 10.5) star with a small semi-major axis of a = 0.0377 ± 0.0005 AU. Ephemeris for the system 
is P = 2.2047299 ± 0.0000040 days, mid-transit time E = 2, 453, 790.2593 ± 0.0010 (BJD). Based 
on the available spectroscopic data on the host star and photometry of the system, the planet has 



a mass of M„ 



1.78 



+ ao5 ^Jup an d radius of R p 



evolved F6 star with M* 



1.47t°;°*M , 



R* = 1.84 



i oc+0.20 n 
1 - ou -0.09 -"-Jup 
+0.23 



-0.11 



oo9-"-Jup- The parent star is a slightly 
i? Qj T off = 6350 ± 80 K, and metallicity 



[Fc/H] = +0.26 ± 0.08. The relatively hot and large host star, combined with the close orbit of the 
planet, yield a very high planetary irradiance of (4.7llJ' 05 ) x 10 9 erg cm~ 2 s _1 , which places the 
planet near the top of the pM class of irradiated planets as defined by Fortney et al. (2007). If as 
predicted by Fortney et al. (2007) the planet re-radiates its absorbed energy before distributing it to 
the night side, the day-side temperature should be about (2730+ 100 ) K. Because the host star is quite 
bright, measurement of the secondary eclipse should be feasible for ground-based telescopes, providing 
a good opportunity to compare the predictions of current hot Jupiter atmospheric models with the 
observations. Moreover, the host star falls in the field of the upcoming Kepler mission; hence extensive 
space-borne follow-up, including not only primary transit and secondary eclipse observations but also 
asteroseismology, will be possible. 

Subject headings: planetary systems — stars: individual (HAT-P-7, GSC 03547-01402) techniques: 
spectroscopic 



1. INTRODUCTION 

Transiting planets provide unique information about 
the nature and evolution of extrasolar planets because 
they yield direct measurement of the radius and (together 
with radial velocity data) mass of these objects. Many 
of these planets are found in quite close-in orbits of their 
parent stars, in which case the radiation from the parent 
star may be expected to play a major role in controlling 
the structure and dynamics of their atmospheres. Tran- 
siting planets with very tight orbits about their parent 
stars are heated sufficiently that their thermal emission 
is measurable, both at times of secondary eclipse and for 
some, throughout their orbit (see e.g. Doming ct al. 2005; 
Knutson et al. 2007; Harrington et al. 2007). Such ob- 
servations, together with observations in different spec- 
tral bands of the depth and shape of primary transits, 
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open up the possibility of detailed understanding of the 
composition, structure, and dynamics of "hot Jupiter" 
atmospheres. (See e.g. Fortney et al. (2007), Hubeny et 
al. (2003), Burrows et al. (2006). 

Here we report on the HATNet project's detection of 
a very hot Jupiter that transits a bright, V = 10.5 F6 
star, with a 2.2 day period and semi-major axis of only 
0.038 AU. The parent star (hereafter denoted HAT-P- 
7), with T e g — 6350 K and radius 1.8 Rq has such a 
high luminosity as to make the surface temperature of 
this very close-orbiting planet among the highest of any 
known transiting planet. 

Because HAT-P-7 is relatively bright, this opens up 
the possibility of a number of interesting follow-on stud- 
ies from space and ground. Notably, the host star is 
also in the field of view of the forthcoming Kepler mis- 
sion. This means that many transits will very likely 
be recorded with extremely high photometric precision 
during the operational phase of Kepler starting in 2009. 
Timing variations in the transits, if detected, could indi- 
cate the presence of other bodies in the system, such as 
terrestrial-mass planets or Trojan satellites. For this rea- 
son, preparatory observations of transits and secondary 
eclipses in the HAT-P-7 system undertaken in advance 
of the mission could significantly increase the scientific 
return from later observations of the star by Kepler. 

The next section (§ 2) describes the photometric detec- 
tion of the HAT-P-7b planet by the HATNet monitoring 
system; § 3 outlines both the photometric and spectro- 
scopic follow-up observations; § 4 explains our analysis of 
these data including derivation of stellar and planetary 
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parameters as well as verification that the signal derives 
from a transiting planet rather than a "false-positive" ; 
and finally, § 5 discusses the results and implications for 
further studies of this intriguing system. 

2. PHOTOMETRIC DETECTION 

The HATNet telescopes HAT- 7 and HAT-8 (HATNet; 
Bakos et al. 2002, 2004) observed HATNet field G154, 
centered at a = 19 h 12 m , S = +45°00', on a near-nightly 
basis from 2004 May 27 to 2004 August 6. Exposures 
of 5 minutes were obtained at a 5.5-minute cadence 
whenever conditions permitted; all in all 5140 expo- 
sures were secured, each yielding photometric measure- 
ments for approximately 33,000 stars in the field down 
to I ~ 13.0. The field was observed in network mode, 
exploiting the longitude separation between HAT-7, sta- 
tioned at the Smithsonian Astrophysical Observatory's 
(SAO) Fred Lawrence Whipple Observatory (FLWO) in 
Arizona (A = 111°W), and HAT-8, installed on the 
rooftop of SAO's Submillimcter Array (SMA) building 
atop Mauna Kea, Hawaii (A = 155° W). We note that 
each light curve obtained by a given instrument was 
shifted to have a median value to be the same as cat- 
alogue magnitude of the appropriate star, allowing to 
merge light curves acquired by different stations and/or 
detectors. 

Following standard frame calibration procedures, as- 
trometry was performed as described in Pal & Bakos 
(2006), and aperture photometry results were subjected 
to External Parameter Decorrelation (EPD, described 
briefly in Bakos et al. 2007), and also to the Trend 
Filtering Algorithm (TFA; Kovacs et al. 2005). We 
searched the light curves of field G154 for box-shaped 
transit signals using the BLS algorithm of Kovacs et al. 
(2002). A very significant periodic dip in brightness was 
detected in the I w 9.85 magnitude star CSC 03547- 
01402 (also known as 2MASS 19285935+4758102; a = 
19 h 28 m 59 s .35, 6 = +47°58'10".2; J2000), with a depth 
of ~ 7.0 mmag, a period off = 2.2047 days and a relative 
duration (first to last contact) of q w 0.078, equivalent 
to a duration of Pq 4.1 hours. 

In addition, the star happened to fall in the overlap- 
ping area between fields G154 and G155. Field G155, 
centered at a = 19 h 48 m , S = +45°00', was also ob- 
served over an extended time in between 2004 July 27 and 
2005 September 20 by the HAT-6 (Arizona) and HAT-9 
(Hawaii) telescopes. We gathered 1220 and 10260 data- 
points, respectively (which independently confirmed the 
transit), yielding a total number of 16620 data-points. 
The combined HATNet light curve after the EPD and 
TFA procedures is plotted on Figs. la,b, superimposed 
on these plots is our best fit model (see § 4) . We note that 
TFA was run in signal reconstruction mode, i.e. system- 
atics were iteratively filtered out from the observed time 
series assuming that the underlying signal is a trapeze- 
shaped transit (see Kovacs et al. 2005, for additional de- 
tails). 

In addition to having a significant overlap with each 
other, we note that fields G154 and G155 both intersect 
the field of view of the upcoming Kepler mission (Borucki 
et al. 2007). 

3. FOLLOW-UP OBSERVATIONS 



3.1. Low-resolution Spectroscopy 

Following the HATNet photometric detection, HAT-P- 
7 (then a transit candidate) was observed spectroscopi- 
cally with the CfA Digital Speedometer (DS, see Latham 
1992) at the FLWO 1.5 m Tillinghast reflector, in order 
to rule out a number of blend scenarios that mimic plan- 
etary transits (e.g. Brown 2003; O 'Donovan et al. 2007), 
as well as to characterize the stellar parameters, such as 
surface gravity, effective temperature, and rotation. Four 
spectra were obtained over an interval of 29 days. These 
observations cover 45 A in a single echelle order centered 
at 5187 A, and have a resolving power of A/AA w 35,000. 
Radial velocities were derived by cross-correlation, and 
have a typical precision of lkms -1 . Using these mea- 
surements, we have ruled out an unblended companion 
of stellar mass (e.g. an M dwarf orbiting an F dwarf), 
since the radial velocities did not show any variation 
within the uncertainties. The mean heliocentric radial 
velocity of HAT-P-7 was measured to be —11 kms -1 . 
Based on an analysis similar to that described in Tor- 
res et al. (2002), the DS spectra indicated that the host 
star is a slightly evolved dwarf with logg = 3.5 (cgs), 
T ff = 6250K and vsini w 6kms _1 . 

3.2. High resolution spectroscopy 

For the characterization of the radial velocity varia- 
tions and for the more precise determination of the stellar 
parameters, we obtained 8 exposures with an iodine cell, 
plus one iodine-free template, using the HIRES instru- 
ment (Vogt et al. 1994) on the Keck I telescope, Hawaii, 
between 2007 August 24 and 2007 September 1. The 
width of the spectrometer slit was 0'.'86 resulting a re- 
solving power of A/AA w 55,000, while the wavelength 
coverage was ~ 3800 — 8000 A. The iodine gas absorp- 
tion cell was used to superimpose a dense forest of I2 lines 
on the stellar spectrum and establish an accurate wave- 
length fiducial (see Marcy & Butler 1992). Relative ra- 
dial velocities in the Solar System barycentric frame were 
derived as described by Butler et al. (1996), incorporat- 
ing full modeling of the spatial and temporal variations 
of the instrumental profile. The final radial velocity data 
and their errors are listed in Table 1. The folded data, 
with our best fit (see § 4.2) superimposed, are plotted in 
Fig. 2a. 

3.3. Photometric follow-up observations 

Partial photometric coverage of a transit event of HAT- 
P-7 was carried out in the Sloan z-band with the Kepler- 
Cam CCD on the 1.2 m telescope at FLWO, on 2007 
November 2. The total number of frames taken from 
HAT-P-7 was 514 with cadence of 28 seconds. During the 
reduction of the KeplerCam data, we used the following 
method. After bias and flat calibration of the images, 
an astrometric transformation (in the form of first or- 
der polynomials) between the ~ 450 brightest stars and 
the 2MASS catalog was derived, as described in Pal & 
Bakos (2006), yielding a residual of ~ 1/4 pixel. Aper- 
ture photometry was then performed using a series of 
apertures of with the radius of 4, 6 and 8 pixels in fixed 
positions calculated from this solution and the actual 
2MASS positions. The instrumental magnitude trans- 
formation was obtained using ~ 350 stars on a frame 
taken near culmination of the field. The transformation 
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TABLE 1 

Relative radial velocity measurements 
of HAT-P-7 
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Fig. 1. — (a) The complete light curve of HAT-P-7 with all of the 
16620 points, unbinncd instrumental /-band photometry obtained 
with four telescopes of HATNet (see text for details), and folded 
with the period of P = 2.2047299 days (the result of a joint fit to 
all available data, § 4.2). The superimposed curve shows the best 
model fit using quadratic limb darkening, (b) Same as a), with 
the transit zoomed-in (3150 data points are shown), (c) Same as 
b), with the points binned with a bin size of 0.004 in phase, (d) 
Unbinncd instrumental Sloan z-band partial transit photometry 
acquired by the KeplerCam at the FLWO 1.2 m telescope on 2007 
November 2; superimposed is the best-fit transit model light curve, 
(e) The difference between the KeplerCam observation and model 
(same vertical scale as in panel d). 



fit was initially weighted by the estimated photon- and 
background-noise error of each star, then the procedure 
was repeated weighting by the inverse variance of the 
light curves. From the set of apertures we have cho- 
sen the aperture for which the out-of-transit (OOT) rms 
of HAT-P-7 was the smallest; the radius of this aper- 
ture is 6 pixels. This light curve was then de-correlated 
against trends using the OOT sections, which yielded a 
light curve with an overall rms of 1.9mmag at a cadence 
of one frame per 28 seconds. This is a bit larger than the 
expected rms of 1.5mmag, derived from the photon noise 
(1.2mmag) and scintillation noise - which has an ex- 
pected amplitude of 0.8mmag, based on the observational 
conditions and the calculations of Young (1967) - pos- 
sibly due to unresolved trends and other noise sources. 
The resulting light curve is shown in Fig. Id. 

4. ANALYSIS 

The analysis of the available data was done in five 
steps. First, an independent analysis was performed 
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on the HATNet, the radial velocity (RV) and the high 
precision photometric follow-up (FU) data, respectively. 
Analysis of the HATNet data yielded an initial value for 
the orbital period and transit epoch. The initial period 
and epoch were used to fold the RV's, and phase them 
with respect to the predicted transit time for a circular 
orbit. The HATNet and the RV epochs together yield a 
more accurate period, since the time difference between 
the discovery light curve and the RV follow-up is fairly 
long; more than 3 years. Using this refined period, we 
can extrapolate to the expected center of the KeplerCam 
partial transit, and therefore obtain a fit for the two re- 
maining key parameters describing the light curve: a/R± 
where a is the semi-major axis for a circular orbit, and 
the impact parameter b = (a/i?*) cos i, where i is the 
inclination of the orbit. 

Second, using as starting points the initial values as de- 
rived above, we performed a joint fit of the HATNet, RV 
and FU data, i.e. fitting all of the parameters simultane- 
ously. The reason for such a joint fit is that the three sep- 
arate data-sets and the fitted parameters are intertwined. 
For example, the epoch (depending partly on the RV fit) 
has a relatively large error, affecting the extrapolation of 
the transit center to the KeplerCam follow-up. 

In all of the above procedures, we used the downhill 
simplex method to search for the best fit values and the 
method of refitting to synthetic data sets to find out 
the error of the adjusted parameters. The latter method 
yields a Monte-Carlo set of the a posteriori distribution 
of the fit parameters. 

The third step of the analysis was the derivation of 
the stellar parameters, based on the spectroscopic anal- 
ysis of the host star (high resolution spectroscopy using 
Keck/HIRES), and the physical modeling of the stellar 
evolution, based on existing isochrone models. As the 
fourth step, we then combined the results of the joint 
fit and stellar parameter determination to determine the 
planetary and orbital parameters of the HAT-P-7b sys- 
tem. Finally, we have done a blend analysis to confirm 
that the orbiting companion is a hot Jupiter. In the fol- 
lowing subsections we give a detailed description of the 
above main steps of the analysis. 

4.1. Independent fits 

For the independent fit procedure, we first analyzed 
the HATNet light curves, as observed by the HAT-6, 
HAT-7, HAT-8 and HAT-9 telescopes. Using the initial 
period and transit length from the BLS analysis, we fit- 
ted a model to the 214 cycles of observations spanned 
by all the HATNet data. Although at this stage we 
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were interested only in the epoch and period, we have 
used the transit light curve model with the assump- 
tion quadratic limb darkening, where the flux decrease 
was calculated using the models provided by Mandel 
& Agol (2002). In principle, fitting the epoch and pe- 
riod as two independent variables is equivalent to fitting 
the time instant of the centers of the first and last ob- 
served individual transits, T c .fi rst and T C; i ast , with a con- 
straint that all intermediate transits are regularly spaced 
with period P. Note that this fit takes into account 
all transits that occurred during the HATNet observa- 
tions, even though it is described only by T c fi rst and 
T clast . The fit yielded T c first = 2453153.0924 ± 0.0021 
(BJD) and T cMst = 2453624.9044 ± 0.0023 (BJD). The 
period derived from the T Cj fi rst and T c j ast epochs was 
PW = 2.20480 ± 0.00049 days. Using these values, we 
found that there were 326 cycles between T Ci i ast and the 
end of the RV campaign. The epoch extrapolated to 
the approximate time of RV measurements was T Ci rv = 
2454343.646 ± 0.008 (BJD). Note that the error in T C;RV 
is much smaller than the period itself (~ 2.2 days), so 
there is no ambiguity in the number of elapsed cycles 
when folding the periodic signal. 

We then analyzed the radial velocity data in the follow- 
ing way. We defined the N tI = transit as that being 
closest to the end of the radial velocity measurements. 
This means that the first transit observed by HATNet 
(at T c , f; rs t) was the iV tri fi rs t = —540 event. Given the 
short period, we assumed that the orbit has been circu- 
larized (Hut 1981) (later verified; see below). The orbital 
fit is linear if we choose the radial velocity zero-point 7 
and the amplitudes A and B as adjusted values, namely: 



2tt 



2tt 



v{t) = 7 + A cos ( — (i-t ) ) +Bsin ( —(*-*„) 

where to is arL arbitrary time instant (chosen to be to = 
2454342.6 BJD), K = \J A 2 + B 2 is the semi-amplitude 
of the RV variations, and P is the initial period P^ 
taken from the previous independent HATNet fit. The 
actual epoch can be derived from the above equation 
since for circular orbits transit center occurs when the RV 
curve has the most negative slope. Using the equations 
above, we derived the initial epoch of the N tT — transit 
center to be T c = 2454343.6462 ± 0.0042 = T^j (BJD). 
We also performed a more general (non-linear) fit to the 
RV in which we let the eccentricity float. This fit yielded 
an eccentricity consistent with zero, namely e cos lu — 
-0.003 ± 0.007 and esinw = 0.000 ± 0.010.' Therefore, 
we adopt a circular orbit in the further analysis. 

Combining the RV epoch T^q with the first epoch ob- 
served by HATNet (T c .fi rst ), we obtained a somewhat re- 
fined period, P^ = 2.204732 ± 0.000016 days. This was 
fed back into phasing the RV data, and we performed 
the RV fit again to the parameters 7, A and B. The fit 
yielded 7 = -37.0± 1.5ms- 1 , K = VA 2 + B 2 = 213.4± 
2.0ms- 1 and = 2454343.6470 ± 0.0042 (BJD). This 

epoch was used to further refine the period to get P^ = 
2.204731 ± 0.000016 d, where the error calculation as- 
sumes that T c q and T c 540 are uncorrelated. At this 
point we stopped the above iterative procedure of refining 
the epoch and period; instead a final refinement of epoch 
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Fig. 2. — (a) Radial-velocity measurements from Keck for HAT- 
P-7, along with an orbital fit, shown as a function of orbital phase, 
using our best fit as period (see § 4.2). The center-of-mass velocity 
has been subtracted, (b) Phased residuals after subtracting the 
orbital fit (also see § 4.2). The rms variation of the residuals is 
about 3.8ms -1 . (c) Bisector spans (BS) for the 8 Keck spectra 
plus the single template spectrum, computed as described in the 
text. The mean value has been subtracted. Due to the relatively 
small errors comparing to the RV amplitude, the vertical scale on 
the (b) and (c) panels differ from the scale used on the top panel. 

and period was obtained through performing a joint fit, 
(as described later in § 4.2). We note that in order to 
get a reduced chi-square value near unity for the radial 
velocity fit, it was necessary to quadratically increase the 
noise component with an amplitude of 3.8 ms , which 
is well within the range of stellar jitter observed for late 
F stars; see Butler et al. (2006). 

Using the improved period P^ and the epoch T c .o, 
we extrapolated to the center of KeplerCam follow-up 
transit (N tv — 29). Since the follow-up observation only 
recorded a partial event (see Fig. Id), this extrapola- 
tion was necessary to improve the light curve modeling. 
For this, we have used a quadratic limb-darkening ap- 
proximation, based on the formalism of Mandel & Agol 
(2002). The limb-darkening coefficients were based on 
the results of the SME analysis (notably, T e g; see § 4.3 

for further details), which yielded 7^ = 0.1329 and 

(z) 

72 = 0.3738. Using these values and the extrapolated 
time of the transit center, we adjusted the light curve 
parameters: the relative radius of the planet p = R p /R±, 
the square of the impact parameter b 2 and the quan- 
tity C/R* = (a/i?*)(27r/P)(l - b 2 )- 1 / 2 as independent 
parameters (see Bakos et al. 2007, for the choice of pa- 
rameters). The result of the fit was p = 0.0762 ± 0.0012, 
b 2 = 0.205±0.144 and (/R+ = 13.60±0.83day-\ where 
the uncertainty of the transit center time due to the rela- 
tively high error in the transit epoch T Cj0 was also taken 
into account in the error estimates. 

4.2. Joint fit 

The results of the individual fits described above pro- 
vide the starting values for a joint fit, i.e. a simultaneous 
fit to all of the available HATNet, radial velocity and 
the partial follow-up light curve data. The adjusted pa- 
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TABLE 2 
Stellar parameters for HAT-P-7 



Parameter Value Source 



T cS (K) 

[Fe/H] 


6350 ± 80 


SME a 




+0.26 ± 0.08 


SME 




v sin i (km s 1 ) . . 


3.8 ±0.5 


SME 




M* (Mq) 


i 47 +0.08 


Y 2 +LC- 


-SME b 


il* (-R©) 


1 84+ ' 23 


Y 2 +LC- 


rSME 


log 9* (cgs) 


4 07 +0.04 
t-vt _ Q8 


Y 2 +LC- 


hSME 


L* (Lq) 


4 9+1.5 


Y 2 +LC- 


rSME 


My (mag) 


3.00 ±0.22 


Y 2 +LC±SME 


Age (Gyr) 


2.2 ± 1.0 


Y 2 +LC- 


-SME 


Distance (pc) 


320+ 50 


Y 2 +LC- 


-SME 



a SME = 'Spectroscopy Made Easy' package for analysis 
of high-resolution spectra Valcnti & Piskunov (1996). Sec 
text. 

b Y 2 +LC+SME = Yalc-Yonsci isochrones (Yi et al. 
2001), light curve parameters, and SME results. 



rameters were T C: _5 40 , the time of first transit center in 
the HATNet campaign, m, the out-of-transit magnitude 
of the HATNet light curve in /-band and the previously 
defined parameters of 7, A, B, p, b 2 and C/R*- We note 
that in this joint fit all of the transits in the HATNet 
light curve have been adjusted simultaneously, tied to- 
gether by the constraint of assuming a strictly periodic 
signal; the shape of all these transits were characterized 
by p, b 2 and (,/ R± (and the limb-darkening coefficients) 
while the distinct transit center time instants were inter- 
polated using T c — 540 — ^c, first and A : B via the KV fit. 
For initial values we used the results of the independent 
fits (§ 4.1). The error estimation based on method refit- 
ting to synthetic data sets gives the distribution of the 
adjusted values, and moreover, this distribution can be 
used directly as an input for a Monte-Carlo parameter de- 
termination for stellar evolution modeling, as described 
later (§ 4.3). 

Final results of the joint fit were: T c -540 = 
2453153.0924 ± 0.0015 (BJD), m = 9.85053 ± 
0.00015mag, 7 = ~37.0il.5ms- 1 , A = 33.8i0.9ms- 1 , 
B = 210.7 i 1.9ms- 1 , p = 0.0763 i 0.0010, b 2 = 
0.135^;^ and C/#* = 13.34 i 0.23 day- 1 . Using the 
distribution of these parameters, it is straightforward to 
obtain the values and the errors of the additional pa- 
rameters derived from the joint derived fit, namely T c ,o, 
a/R±, K and P. All final fit parameters are listed in 
Table 3. 

4.3. Stellar parameters 

The results of the joint fit enable us to refine the pa- 
rameters of the star. First, the iodine-free template spec- 
trum from Keck was used for an initial determination of 
the atmospheric parameters. Spectral synthesis model- 
ing was carried out using the SME software (Valcnti & 
Piskunov 1996), with wavelength ranges and atomic line 
data as described by Valcnti & Fischer (2005). We ob- 
tained the following initial values: effective temperature 
6350 i 80 K, surface gravity log g* = 4.06 i 0.10 (cgs), 
iron abundance [Fc/H] = +0.26 i 0.08, and projected 
rotational velocity vsini — 3.8 i 0.5kms _1 . The rota- 
tional velocity is slightly smaller than the value given 
by the DS measurements. The temperature and surface 
gravity correspond to a slightly evolved F6 star. The 



uncertainties quoted here and in the remaining of this 
discussion are twice the statistical uncertainties for the 
values given by the SME analysis. This reflects our at- 
tempt, based on prior experience, to incorporate system- 
atic errors (e.g. Noyes et al. (2008); see also Valenti & 
Fischer (2005)). Note that the previously discussed limb 

darkening coefficients, "f[ z \ 72 Z \ an d 7^ have been 
taken from the tables of Claret (2004) by interpolation 
to the above-mentioned SME values for T c ff, logg*, and 
[Fe/H] . 

As described by Sozzetti et al. (2007), a/i?* is a bet- 
ter luminosity indicator than the spectroscopic value of 
log g+ since the variation of stellar surface gravity has 
a subtle effect on the line profiles. Therefore, we used 
the values of T e s and [Fe/H] from the initial SME anal- 
ysis, together with the distribution of a/i?* to estimate 
the stellar properties from comparison with the Yonsei- 
Yale (Y 2 ) stellar evolution models by Yi et al. (2001). 
Since a Monte-Carlo set for a/R± values has been de- 
rived during the joint fit, we performed the stellar pa- 
rameter determination as follows. For a selected value 
of a/R±, two Gaussian random values were drawn for 
T e g and [Fe/H] with the mean and standard deviation as 
given by SME (with formal SME uncertainties doubled 
as indicated above). Using these three values, we searched 
the nearest isochrone and the corresponding mass by us- 
ing the interpolator provided by Demarque et al. (2004). 
Repeating this procedure for values of a/R±, T ff, [Fe/H], 
the set of the a posteriori distribution of the stellar pa- 
rameters was obtained, including the mass, radius, age, 
luminosity and color (in multiple bands). The age deter- 
mined in this way is 2.2 Gy with a statistical uncertainty 
of i0.3 Gy; however, the uncertainty in the theoretical 
isochrone ages is about 1.0 Gy. Since the corresponding 
value for the surface gravity of the star, log 5* = 4.07l , ' '| 
(cgs), is well within l-er of the value determined by the 
SME analysis, we accept the values from the joint fit as 
the final stellar parameters. These parameters are sum- 
marized in Table 2. 

We note that the Yonsei-Yale isochrones contain the 
absolute magnitudes and colors for different photometric 
bands from U up to M , providing an easy comparison 
of the estimated and the observed colors. Using these 
data, we determined the V — I and J — K colors of the 
best fitted stellar model: (V - I) YY = 0.54 i 0.02 and 
(J - K) YY = 0.27 i 0.02. Since the colors for the in- 
frared bands provided by Yi et al. (2001) and Demarque 
et al. (2004) are given in the ESO photometric standard 
system, for the comparison with catalog data, we con- 
verted the infrared color (J — K)yy to the 2MASS sys- 
tem (J — Ks) using the transformations given by Car- 
penter (2001). The color of the best fit stellar model 
was (J — Ks)yy = 0.25 i 0.03, which is in fairly good 
agreement with the actual 2MASS color of HAT-P-7: 
(J - Ks) = 0.22 i 0.04. We have also compared the 
(V — I)yy color of the best fit model to the catalog data, 
and found that although HAT-P-7 has a low galactic lat- 
itude, 611 = 13?8, the model color agrees well with the 
observed TASS color of (V - I)tass = 0.60 i 0.07 (see 
Droege et al. 2006). Hence, the star is not affected by the 
interstellar reddening within the errors, since E{V — I) = 
(V - I)tass - (V- I)yy = 0.06 i 0.07. For estimating 
the distance of HAT-P-7, we used the absolute magnitude 
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My = 3.00 ± 0.22 (resulting from the isochrone analy- 
sis, see also Table 2) and the Vtass = 10.51 ± 0.06 ob- 
served magnitude. These two yield a distance modulus of 
Vtass -My = 7.51 ±0.28, i.e. distance of d = 320±^pc. 

4.4. Planetary and orbital parameters 

The determination of the stellar properties was fol- 
lowed by the characterization of the planet itself. Since 
Monte-Carlo distributions were derived for both the light 
curve and the stellar parameters, the final planetary and 
orbital data were also obtained by the statistical analysis 
of the a posteriori distribution of the appropriate combi- 
nation of these two Monte-Carlo data sets. We found 



TABLE 3 
Orbital and planetary parameters 



that the mass of the planet is M p = 1.776 



the radius is R v 



1.363 



+0.195 



+0.077 

-0.049 



M, 



-0.087 an d its density is 
p p = 0.876^^24 gcm~ 3 . We note that in the case of 
binary systems with large mass and radius ratios (such 
as the one here) there is a strong correlation between 
M p and R p (see e.g. Beatty et al. 2007). This correlation 
is also exhibited here with C(M p , R p ) — 0.89. The final 
planetary parameters are also summarized at the bottom 
of Table 3. 

Due to the way we derived the period, i.e. P = 
(T c ,o — T C! _54o)/540, one can expect a large correlation 
between the epochs T c o, T c -540 and the period itself. 
Indeed, C(T C _ 540 ,P) = -0.783 and C(T cfi ,P) = 0.704, 
while the correlation between the two epochs is negligi- 
ble; C(T C _54o,T Ci0 ) = -0.111. It is easy to show that 
if the signs of the correlations between two epochs Ta 
and Tb (in our case T c , and T C! _ 5 4o) and the period 
are different, respectively, then there exists an optimal 
epoch E, which has the smallest error among all of the 
interpolated epochs. We note that E will be such that it 
also exhibits the smallest correlation with the period. If 
ct(Ta) and <t(Tb) are the respective uncorrelated errors 
of the two epochs, then 

r r A q(TB) 2 + r B q(r A ) 2 l 

[ ^Te^^Ta) 2 \ [) 

where square brackets denote the time of the transit 
event nearest to the time instance t. In the case of 
HAT-P-7b, T A = T c _ 540 and T B = T Cj0 , the correspond- 
ing epoch is the event N tv = —251 at E = T c ,_25i = 
2,453, 790.2593 ±0.0010 (BJD). The final ephemeris and 
planetary parameters are summarized in Table 3. 

4.5. Excluding blend scenarios 

Following Torres et al. (2007), we explored the possi- 
bility that the measured radial velocities are not real, but 
instead caused by distortions in the spectral line profiles 
due to contamination from a nearby unresolved eclipsing 
binary. In that case the "bisector span" of the aver- 
age spectral line should vary periodically with amplitude 
and phase similar to the measured velocities themselves 
(Queloz et al. 2001; Mandushev et al. 2005). We cross- 
correlated each Keck spectrum against a synthetic tem- 
plate matching the properties of the star (i.e. based on 
the SME results, see § 4.3), and averaged the correlation 
functions over all orders blueward of the region affected 
by the iodine lines. From this representation of the aver- 
age spectral line profile we computed the mean bisectors, 
and as a measure of the line asymmetry we computed the 
"bisector spans" as the velocity difference between points 



Parameter 



Value 



Light curve parameters 

P (days) 2.2047299 ± 0.0000040 

E (BJD — 2,400,000) 53, 790.2593 ± 0.0010 

T14 (days) a 0.1685 ± 0.0110 

T12 = T 34 (days) a 0.0150 ± 0.0036 

" 4-3518:11 

Rp/R* 0.0763 ± 0.0010 

6 = acosi/i?* O.37+029 

i (deg) 85?7t|;5 

Spectroscopic parameters 

K (ins" 1 ) 213.5 ± 1.9 

7 (kms- 1 ) -37.0 + 1.5 

e (adopted) 

Planetary parameters 

M P (Afj) 1.776l8;°S 

-Rp (-Rj) 1. 363^0 087 

C(M P ,R P ) 0.89 

Pp (gem" 3 ) 0.8761°-^ 

a (AU) 0.0377 ± 0.0005 

\ogg p (cgs) 3.31 + 0.08 

a Tl 4 : total transit duration, time between first to last 
contact; T\2 = T34: ingress/egress time, time between first 
and second, or third and fourth contact. 

selected near the top and bottom of the mean bisectors 
(Torres et al. 2005). If the velocities were the result of 
a blend with an eclipsing binary, we would expect the 
line bisectors to vary in phase with the photometric pe- 
riod with an amplitude similar to that of the velocities. 
Instead, we detect no variation in excess of the measure- 
ment uncertainties (see Fig. 2c) . We have also tested the 
significance of the correlation between the radial veloc- 
ity and the bisector variations. Therefore, we conclude 
that the velocity variations are real and that the star is 
orbited by a Jovian planet. We note here that the mean 
bisector span ratio relative to the radial velocity ampli- 
tude is the smallest (~ 0.026) among all the HATNet 
planets, indicating an exceptionally high confidence that 
the RV signal is not due to a blend with an eclipsing 
binary companion. 

5. DISCUSSION 

With orbital period of only 2.2 days, implying a semi- 
major axis of only 0.038 AU, HAT-P-7b, is a "very hot 
Jupiter" , a name frequently given to giant planets with 
orbital period less than 3 days (see Bouchy et al. 2004). 
However, what really matters is the incident flux from 
the parent star, which for HAT-P-7b is exceptionally high 
because (a) the large radius of the parent star (1.84 Rq) 
combined with the small semi-major axis yields a very 
small geometrical factor a/R* — 4.35 (whose inverse de- 
scribes the stellar flux impinging on the planet); and 
(b) its parent star has relatively high effective tempera- 
ture: T ff = 6350 K. The flux impinging on the planet is 
4.7 x 10 9 erg cm _2 s _1 . This places it among the most 
highly irradiated members of the "pM" class of planets, 
as discussed by Fortney et al. (2007). 

For a pM class planet, Fortney et al. (2007) argue 
that the incident radiation is absorbed by TiO and VO 
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molecules in a hot stratosphere, and almost immediately 
re-radiated, rather than being partially advected to the 
night side of the planet, as would be the case if the radia- 
tion were absorbed deeper down in the cooler atmosphere 
of a "pL" planet, which would not have gaseous TiO and 
VO in its atmosphere. 

For a pM planet, the dayside temperature would be 
given by 

T p = T eff (^y /2 [/(l-^)] 1/4 - (3) 

where Ab is the Bond albedo, which we assume to be 
essentially zero (see, e.g. Rowe et al. 2007), and the re- 
distribution factor / = 2/3 (see Lopez-Morales and Sea- 
ger 2007). For HAT-P-7 this yields a dayside tempera- 
ture of (2730^}qq) K. If at the other extreme, contrary to 
the theory of Fortney et al. (2007), the incident energy 
absorbed by HAT-P-7b were to be distributed isotrop- 
ically over the planet's surface before being re-emitted, 
as is predicted for "pL" planets too cool to have signif- 
icant TiO and VO in their atmospheres, then / = 1/4 
in Equation 3, and the day-side temperature would be 
(2140^eo°) K. Thus, detection and measurement of the 
emission from the dayside of HAT-P-7b would yield an 
important test of the theory of Fortney et al. (2007) . 

Lopez-Morales and Seager (2007) and Fortney et al. 
(2007) noted that the thermal emission in the optical and 
near infrared from very hot Jupiters should be detectable 
from the ground. Lopez-Morales and Seager (2007) iden- 
tified OGLE-TR-56b, the only other known very hot 
jupiter with predicted dayside temperature as high as 
that of HAT-P-7b, as a particularly promising candi- 
dates. However, since HAT-P-7 is 6 magnitudes brighter 
than OGLE-TR-56, it is a far better candidate for sec- 
ondary eclipse studies from the ground. For HAT-P-7, 
the brightness decrease in z'-band at secondary eclipse 
should be ~ 0.2 mmag or ~ 0.04 mmag in the respective 
cases of zero (pM-typc) or complete (pL-type) planet- 
wide re-distribution of irradiated stellar flux. This im- 
plies that in the former case the eclipse could be detected 
at the 2-sigma level with a 1.2-m telescope with a z'-band 
integration of 4 hours during eclipse (that is, the total 
eclipse duration), plus 4 more hours of out-of-eclipse ob- 
servation. We plan to carry out such observations in the 
near future with the KeplerCam on FLWO 1.2 m and 
with MMT 6.5 m (also at FLWO) telescopes. 

Recently Pollacco et al. (2007) reported the discovery 
of the transiting exoplanet WASP-3b, located 0.0317 AU 
from an F star with T g = 6350 K and radius i?* = 
1.31i?0; this indicates a very high radiative flux at the 
planet of 2.5 x 10 9 erg cm~ 2 s -1 , comparable to but about 
25% smaller than the radiative flux at HAT-P-7b. For 
the same reasons as for the case of HAT-P-7b discussed 
above, this almost certainly will be a pM-class planet, 
so that its dayside temperature should be about 2530K, 
as indicated by Equation 3 with / = 2/3 and Ab = 0. 
This planet should also be detectable in visible or near 
infrared radiation. In fact, if it is a pM planet, the depth 
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of secondary eclipse as seen in z'-hand should be similar 
to that for HAT-P-7b, namely ~ 0.2 mmag: although 
with its lower dayside temperature the surface brightness 
of WASP-3b should be smaller than for HAT-P-7b, the 
geometrical factor (i? p /i?*) 2 is larger, because WASP-3 
is a less evolved star. 

HAT-P-7b has an additional important property in 
that it is in the field of view of one of the detectors on 
the forthcoming Kepler mission (Borucki et al. 2007), 
currently scheduled for launch in early 2009 to moni- 
tor stars for photometric signals from transiting planets 
and other forms of stellar variability. Hence, HAT-P-7 
should be subjected to intense and regular photometric 
observations for several years starting in 2009. This will 
permit extremely high precision light curves of the pri- 
mary transit, the secondary transit, and perhaps also the 
variation of thermal and reflected light emission from the 
planet over the course of its orbit as described above. In 
addition, it will be possible to search intensively for tran- 
sit timing variations over the several year lifetime of the 
mission. Such observations could yield important infor- 
mation about orbital variations due to other companions 
in the system including terrestrial-mass or smaller plan- 
ets or Trojans, or to other effects such as precession of 
the orbital plane or the orbital line of nodes. To maxi- 
mize the scientific return from Kepler in this regard, it 
would be very useful to obtain high precision light curves 
both from the ground, and from space instruments such 
as Epoxi and MOST, in the very near term for later com- 
parison with Kepler observations. 

Asteroseismology of HAT-P-7 using Kepler could yield 
an independent measurement of the star's mean density 
from the "large separation" Av between p-modes of the 
same low angular degree I (e.g. Stello et al. 2007). This 
would provide a check on the mean density of the star 
completely independent of that derived from the transit 
measurements of the a/ R± parameter. Equivalently, the 
asteroseismology measurements would yield an indepen- 
dent measurement of the stellar radius, and hence a more 
precise value for the radius of the planet HAT-P-7b. Such 
a rich combined data set could have additional benefits 
for understanding stellar structure and evolution. 
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